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(S7> Abstract- A device for noninvasive, continuous moniloring of arterial blood pressure for advanced cardiovascular «ii»gn^ 
St SSt^t^^ive. continuom blood p«ssure n^easuxement devices are mechanically intrusive and, 

f«tn™ ambulatory monitoring. Ibis new approach mquires only simple, noninvasive '-"j-^^ de^-^s«^^fi^^ 
jSorfeth5^ographs and anelectrical impedance photoplethysmograph (EIP) to monitor the dynamic behavior of the ar^^ Wood 
floS to S^Si. .neasmed signalsfi^ttese noninvasive sensors on an arterial segment are mtegrated oestuna^the blood 
o^;m^S^StSl on a hemodynamic model. A mathematical model of the arterial blood flow ,s denved and transfonned 
K a^S-SS^^«^ ^ *e mileling. a p«cise hemodynamic model for the arterial segment on which sensors are loca^ 
S^Swith .datively simplified models of the upstream and the downstr«un artenal flows to represent an entire 
i^SS^ ^TlCahnan filter based «>n tbe model and it is shown that the internal vanables such as die artena^ 

S t>S3' in iaial segment can be estimated based On the measurements, even though the observabdity condition of the 
":^m'^;^t^^r?iSllts indicate th^^ 
in real-time even from noisy sensor signals. 
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CUFFLESS CONTINUOUS BLOOD PRESSURE MONITOR 

t 

* 

Technical Field 

The present invention relates to a device and method for monitoring the blood pressure of 
a patient and, more panicularly, for deriving the blood pressure from measurements performed 
continuously on the finger of the patient. 

Background of the Invention 

« 

Noninvasive ambulatory blood pressure monitoring is currenUy limited to the simple 
measurements of systolic and diastolic blood pressures at intervals. However, it is known to 
clinicians that continuous waveforms of the blood pressure can provide more useful information 
about the patient's cardiovascular state that are difficult to obtain from the routine antecubital 
pressure measurement. For example, the rate of pressure rise at the beginning of systole indicates 
the strength of cardiac contraction while the rate of pressure decay during end diastole can be 
used as a measure of peripheral vascular resistance, both of which are important parameters used 
in cardiovascular diagnoses. In fact, many numerical algorithms have been developed to estimate 
left-ventricular and circulatory parameters from the arterial pressure waveform by applying a 
computer model of the cardiovascular system, as described by J. W. Clark, et al., "A Two-Stage 
Identification Scheme for the Determination of the Parameters of a Model of the Left Heart and 
Systemic Circulation," IEEE Trans, on Biomed. Eng., Vol. 27. pp. 20-29. Jan.. 1980; W. 
Welkowitz. Q. Cui. Y. Qi and J. Kostis, "Noninvasive Estimation of Cardiac Output," IEEE 
Trans, on Biomed. Eng.. Vol, 38. pp. 1100-1 105, Nov.. 1991; M. Guarini, J. Urzua, A. Cipiiano. 
and W. Gonzalez. "Estimation of Cardiac Funcuon From Computer Analysis of the Arterial 
Pressure Wavefonn," IEEE Trans, on Biomed. Eng.. Vol. 45. pp. 1420-1428. Dec. 1998; and E. 
T. Ozawa. "A Numerical Model of the Cardiovascular System for Clinical Assessment of the 
Hemodynamic State," Ph.D. Thesis, Dept. of Health Sciences and Technology, MIT, Sep.. 1996. 
Considering that heart disease is a prevalent cause of death in the modem society, it is obvious 
that long-term noninvasive continuous monitoring of such pressure waveforms would bring 
enormous improvement of the quality of healthcare at home as well as in the hospital. 

A few devices have been developed for continuous monitoring of the arterial pressure 
wavefonn. yet these are either invasive or mechanically intrasive and are not designed for the 
long-tenn use. For example. Pressman and Newgard developed a noninvasive method for 
continuously measuring the instantaneous blood pressure by applying the coplanar measurement 
principle used by tonometry, as described in G. Pressman and P. Newgard, "A Transducer for 
Conunuous External Measurement of Arterial Blood Pressure." IEEE Trans, on Biomed. Eng.. 
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Vol. 10, pp. 73-81, 1961. In this method, called "arterial tonometry." the artery is flattened by 
applying external pressure non-invasively to squeeze the anery against the bone. Since the 
circumferential tension of the arterial wall disappears, the applied pressure to maintain the 
flattened shape indicates the arterial blood pressure. An array of piezoelectric transducers is used 
5 for the pressure reading. Penaz, on the other hand, proposed a new noninvasive, continuous 
blood pressure measuring method based on the principle of vascular wall unloading, as described 
in Penaz, "Photo-electric Measurement of Blood Pressure, Volume and Bow in the Finger," 
Digest of the 10-th Int. Conf. on Medical and Biolog. Eng., 1973. In this method, a cuff is 
inflated to a pressure equal to the pressure in the artery and the cuff pressure is conunuously 
10 adjusted by a servo control system, which monitors the size of the artery using a 

photoplethysmograph. This meth9d was further developed by Wesseling, as described in K. H. 
Wcsseling, "Non-invasive, Continuous, Calibrated Blood Pressure by the Method of Penaz," 
Blood Pressure Measurement and Systemic Hypertension, pp.163-175. Medical Worid Press, 
and successfully commercialized as "FINAPRES." Yamakoshi and his group also developed a 
15 similar device independenUy by applying the vascular unloading technique, as described in C. 
Tase and A. Okuaki. "Noninvasive Continuous Blood Pressure Measurement - Clinical 
Application of FINAPRES - Japanese J. of Clinical Monitor, Vol. 1, pp.61 -68, 1990; and K. 
Yamakoshi, H. Shimazu and T. Togawa, "Indirect Measurement of Instantaneous Arterial Blood 
Pressure in the Human Finger by the Vascular Unloading Technique," IEEE Trans, on Biomed. 
20 Eng., Vol. 27, pp. 150-155, 1980. The major drawback of these devices, however, is the tight 
confinement and mechanical inuiisiveness of the sensor probes and the resultant discomfort to 
the patient. As stated above, these methods require a constant and conUnuous external pressure 
on the skin surface of the patient and it could cause vasospasm and pressure drops in the 
peripheral artery, as described in A. Kawarada, H. Shimazu, H. Ito. and K. Yamakoshi. 
25 "Ambulatory Monitoring of Indirect Beat-to-Beat Arterial Pressure in Human Fingers by a 
Volume-Compensation Method," Med Biol Eng Comput, Vol. 34, pp. 55-62. Jan. I99I. For 
long-term, ambulatory blood pressure monitoring, a new method for the noninvasive and non- 
intrusive continuous measurement is preferable. 

Summary of the Invention 
30 In accordance with one aspect of the invention, a system for monitoring a blood pressure 

of a patient uses a first phptoplethysmograi)h proximate to a finger of the patient for providing a 
measure of change in the arterial diameter at a first location of a specified 'artery of the patient. A 
second photoplethysmograph proximate to the finger of the patient and displaced relative to the 
first photoplethysmograph provides a measure of change in the arterial diameter at a second 
35 location of the specified artery of the patient. An elecuical impedance plethysmograph in 
electrical contact with the finger of the patient provides a measure of change in the electrical 
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impedance of an arterial segment between the first and the second locations of the specified 
artery. A controller derives a measure of the blood pressure of the patient based on the measures 
of change in the arterial diameter at the first and second locations of the specified artery and the 
measure of change in the electrical impedance of an arterial segment. 
5 In a further related embodiment, the first photoplethysmograph is borne by the patient on 

a finger ring. In another related embodiment, the first photoplethysmograph is borne by the 
patient on a first band of a finger ring and the second photoplethysmograph is borne by the 
patient on a second band of the finger ring. In some embodiments, a transmitter may optionally 
be used for transmitting the measure of the blood pressure of the patient to a remote location. 
10 In another embodiment, a system for monitoring a blood pressure of a patient uses a 

monitor having a first and a second band to be worn by the patient on a single finger. The 
monitor has a first photoplethysmograph disposed on the first band for providing a first signal 
based on a first arterial diameter of the patient, a second photoplethysmograph disposed on the 
second band for providing a second signal based on a second arterial diameter of the patient, and 
15' an electrical impedance plethysmograph disposed on the first and second bands for providing a 
third signal based on the electrical impedance of the a segment of an artery of the patient. A 
conuoller analyzes the first, second , and third signals and determines a measure of the blood 

pressure of the patient. 

In accordance with another embodiment, a method for monitoring the blood pressure of a 
20 patient derives a measure of change in both the diameter of the first and second end of a segment 
of an artery of the patient. A model of arterial blood flow is applied to the derived measures of 
change in the diameters of the first and second ends of the arterial segment and the volume of the 
segment for calculating the instantaneous blood pressure of the patient. In a related embodiment, 
the step of deriving the measure of change in the diameter of the first end of a segment of an 
25 artery includes receiving a signal of a first photoplethysmograph. In another related 

embodiment, the step of deriving the measure of change in the diameter of the second end of a 
segment of an artery includes receiving a signal of a second photoplethysmograph. In some 
embodiments, a Kaman filter is used for estimating internal state variables. 

Brief Description of the Drawings 
30 The foregoing features of the invention will be more readily understood by reference to 

the following detailed description, taken with reference to the accompanying drawings, in which: 
Figure 1 schematically shows a Kalman filter for instantaneous blood pressure estimation 
in accordance with an embodiment of the present invention. 

Figure 2 shows a segment of a viscoelastic artery with length of L. 
35 Figure 3 shows a state of stress in a thin-walled, viscoelastic blood vessel. 

Figure 4 shows discretization of the hemodynamic model of a digital arterial sebment. 
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Figure 5 shows an extended Windkessel model for upstream dynamics. 
Figure 6 shows a classic Windkessel model for downstream dynamics. 
Figure 7 shows cuff-less ambulatory pressure monitoring in accordance with an 
embodiment of the present invention. 
5 Figure 8 shows an artery model used for simulation. 

Figure 9 shows a plot of system input: blood pressure on the boundary as a function of 

time. 

Figure 1 Oa shows a plot of system outputs: arterial section areas SI and S3 as a funcuon 
of time. 

10 Figure 1 Ob shows a plot of system outputs: volumetric change V as a ftjncuon of time. 

Figure 1 1 shows a plot of output measurement V and estimation as a function of time. 
Figure 1 2 shows a plot of digital blood pressure estimation by a Kalman Filter vs. 
measurement by an arterial tonometer as a function of time. 
Figure 13 shows an exemplary monitoring system. 
15 Description of Specifi r F.mhodimentS 

A new approach to noninvasive non-introsive continuous measurement of pulsating 
arterial blood pressure is now described which does not require the use of a cuff. In accordance 
with embodimems of the invention, information gathered by multiple sources and sensors are 
merged to provide improved insight into the phenomena under consideration. The use of sensor 
20 fusion is applied for indirecUy estimating arterial blood pressure by integrating simultaneous 
measurements from noninvasive, non-intrusive sensors such as a photoplethysmograph and a 
bioelecuical impedance plethysmograph with a mathematical model of the blood flow. 

The present invention may be applied in the context of finger-sensors such as those 
described in U.S. Pat. No. 5.964,701 which is herein incorporated by reference. 
25 In this approach, a Kalman filter is used for the sensor fusion scheme. Kalman filters arc 

the standard state estimators or observers that are optimum with respect to the process noise and 
sensor noise, and many nonlinear extensions have been developed, and are well known to people 
skilled in the art. as described in R.G. Brown and P. Y.Z. Hwang. Introduction to Random 
Signals and Applied Kalman Filtering, John Wiley and Son. 1997. hereby incorporated by 
30 reference. In accordance with preferred embodiments, a state-space equation is derived from a 
mathematical hemodynamic model and a Kalman filter is applied to esumate the internal stale 
variables such as the blood pressure based on signals from noninvasive and non-intrasive 
sensors. Figure 1 shows the basic scheme of this approach. 

Furthermore, in accordance with preferred embodiments, a two-dimensional 
35 mathematical model of the arterial blood flow is derived as an incompressible, axially symmetric 
Newtonian fluid in a rectilinear, viscoelastic thick shell of isotropic, incompressible material 
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with a circular section. The modeling meihod is applied to a small digital arterial segment, from 
which sensor signals such as a photoplethysmograph and a bioelectrical impedance 
plethysmograph are obtained. Then, the hemodynamic model of the peripheral arterial segment is ' 
extended up to the heart as the proximal boundary and the capillary as the distal boundary to 
5 represent an entire arterial stream. A commonly assumed pattern of the cardiac output is used as 
the system's input. To avoid high-order modeling, the upstream is modeled as a three- 
dimensional Windkessel model, and the downstream is modeled as simple impedance. Finally, a 
Kalman filter is designed based on the extended model. Since the original local arterial segment 
are precisely modeled and the output signals are measured from the segment, it is expected that 
10 the Kalman filter can estimate the local arterial blood pressure accurately even with the 
simplifications of the input and the modeling of the upstream and downstream blood flows. 

Slate-Space Modeling of Arterial Hemodynamics 

15 A mathemaucal hemodynamic model is used for the complex behavior of the arterial vessel 
and blood flow of a peripheral arterial segment and internal variables such as the blood pressure 
are estimated by comparing the sensor readings from the segment with the simulated outputs. 
Therefore, the accuracy and fidelity of the local model is a key issue. Many hemodynamic 
models have been developed for the study of the two-dimensional nonlinear behavior of the 
20 pulsating blood flow, as described in J. C. SietUer, P. Niederer and M. Anliker, "Theoretical 
Analysis of Arterial Hemodynamics including the Influence of Bifurcations," Annals of Biomed. 
Eng., Vol, 9, pp. 145-164, 1981; and G. A. Johnson, H. S. Borovetz, and J. L. Anderson, "A 
Model of Pulsatile Flow in a Uniform Deformable Vessel," J. of Biomechanics, Vol. 25, pp. 91- 
100, 1992. A mathematical framework developed by Belardinelli and Cavalcanti, as described in 
25 E. Belardinelli and S. CavalcanU. "A New Nonlinear Two-Dimensional Model of Blood Motion 
in Tapered and Elastic Vessels," Comput. Biol. Med., Vol. 21, pp. 1-13, 1991; and E. 
Belardinelli and S. Cavalcanti, "Theoretical Analysis of Pressure Pulse Propagation in Arterial 
Vessels," J. of Biomechanics, Vol. 25, pp. 1337-1349, 1992, is applied which describes a two- 
dimensional nonlinear flow of Newtonian viscous fluid moving in a defonnable tapered tube. 
30 The papers of Belardinelli and Cavalcanti arc incorporated herein by reference. The upstream 
and the downstream arterial flows are represented as an extended Windkessel model, and 
combined with the above nonlinear model of the local segment to constitute the entire 
stream. 

35 Local Arterial Flow Model 
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Mathematical Model of Arterial Flow 

A small segment (disiance of L) of a small anery such as a digital artery is shown in 
Figure 2. The arterial vessel is assumed to be a rectilinear, deformable, thick shell of isotropic. . 
incompressible material with a circular section and without longitudinal niovemenis. Blood is an 
5 incompressible Newtonian fluid and flow is axially symmetric. Two-dimensional Navier-Stokes 
equations and continuity equation for a Newtonian and incompressible fluid m cyhndncal 
coordinate (r, ft z) arc: 

du du^ du_ 1 . „.i!ii+l3fi+£fi) 0) 

dw 1 .^(^+L^+£21-2!L^ (2) 

r dr oz 

where P denotes pressure, p density, vkinematic viscosity, and u^u(r,z.t) and w=wfr.z./> denote 
the components of velocity in axial (z) and radial (r) directions respectively, as shown m Fig. 2. 
Let R(z,t) denote the inner radius of the vessel and define a new variable: 

The pressure P is assumed to be uniform within the cross section so that P is independent of the 
radial coordinate. 7J, i.e. />=P(z.O- The above equations can be rewritten in a new coordmate (i?. 



15 



ir"i^^''^"ar^ar^""^a^^"ar-:^%^%a;7 

1 aw w ^ au a/? 9" _ Q (?) 

Rdrj rjR dz Rdzdfj 

where it can be assumed: 

The boundary conditions for the above equations in rj axis are: 
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The basic idea of this hemodynamic modeling, described by E. Belardinelli and S. Cavalcanti. is 
to assume thai the velocity profile in the axial direction can be expressed as the following 
polynomial form: 



The velocity profile in the radial direction is also expressed as: 



10 

For simplicity, N= 1 , such as 

15 

By plugging eqs.(l 1) and (12) into eqs.(5) and (7), the dynamic equations of qUJ) and R{z,i) arc 
obtained as: 

^_l£^_^^ + iL^+li^ = 0 (13) 
dt R dt R dz R^ pdz 

^„dR R^ dg ^ dR ^ ,j4^ 

dt 2 dz dz 

20 

Complete derivaiions of the above equations are described by described by E. BelardinelU and S. 
Cavalcanii. The cross-seciional area 5(2.0 and blood flow Qiz.t) can be defined as: 



s ^ 



25 Then, eqs.( 1 3) and ( 1 4) can be re- written in leims of Q and 5 as: 



+ + — — = 0 (15) 

3/ S dt 5' dz S Ipdz 

^+^=0 (16) 
dr dz 
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Viscoelastic Model of Arterial Wall 

« 

To undersiand ihe hemodynamics of anerial blood flow, a modeling of the viscoelastic behavior 
5 of the arterial wall is essential. A constitutive law of the arterial wall is derived from the stress- 
strain relationship of the material w. Let a^and a, be the circumferential stress and tangential 
stress respectively as shown in Figure 3. Ignoring the inertia of the arterial wall and the external 
pressure, equilibrium with the blood pressure gives: 



PR = a ge - a ,eR 



(17) 



10 



where /?(z,0 and e are the radius of the anerial vessel and the thickness of the arterial wall 
respectively. 

I 

From the geometric compaUbility of the blood vessel, an expression of strains can be obtained 
IS stich as 



ee = 



R-R, 



(18) 



where ft and e are circumferemial and tangential strains respectively and a constant Ro is the 
20 radius of the artery when P(z.r)=0 and the system is in a steady state. 

The most widely used model to describe the viscoelastic properties of the arterial wall is the 
Kelvin-Voigt model, in which the stress-suain relationship is described as: 



25 



O0 = E£0 +7?-—-. cr, =E£,+7J 



dt 



dt 



(19) 



in which E is the elastic modulus and rj is the damping coefficient. By plugging eqs.(18) and 
(19) with S&=JiRo' and eliminaUng second and higher order tenns, the following equation is 
obtained describing the viscoelastic constitutive law of the anerial wall: 



30 



(20) 
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Discretization 

The above nonlinear, partial differential equations given in eqs.(15), (16) and (20) are discretized 
and transformed into a state equation using a finite-difference method. First, the segment of the 
anery (length L) is equally divided by N grids with a step size of Az=U(N-l). The mesh points in 
the finite difference grids are represented by j where j=l,2,- •,N and N > 2. If the length of the 
arterial element Az is sufficienUy small then it is possible to approximate - in each section - the 
derivatives with respect to the axial coordinate z with the following finite difference scheme: 

85,. ^ 5,„ - S, dP, _ Pi.i-Pi ^ 96. ^ g, - e.-. (21) 
dz Az ' dz Az ' dz Az 



The constitutive law given in eq. (20) is modeled such that the viscoelasticity applies only at 
mesh points. An example of the discretization when N=4 is shown in Figure 4. Using the above 
15 equations, the hemodynamic model given in eqs.(15) and (16) can be discreuzed as 

dQ, 3Q, q 2Q,^ 5,., - 5, ^ Axv ^ 5^i^>,-/^ (22) 

dt 5, Az S," Az 5, • 2p Az 

dSj _ Qi - g.-i (23) 
dt Az 

To complete the discretization of the hemodynamic model, the boundary condition at proximal 
20 (P/. Qo) and distal {Pn. Qn) extremities of the arterial segment must be defined appropriately. 

Upstream Blood Flow 

Upstream dynamics extends the proximal boundaiy (Pj, Qo) up to the heart so that the commonly 
25 assumed pattern of the cardiac output can be used as the input to the system. For simplicity, a 
lumped model is used to describe the upstream dynamics- A large amount of work has been done 
in this area. In accordance with a preferred embodiment, a four-element modified Windkessel 
model is applied, as described in G. Landes. Einige uniersuchungen an elektrischen analogie- 
schaltungen zum kreislauf-system. Z. Biol., 101:410, (1943) This model has been adopted by 
30 many researchers for the arterial pressure waveform analysis, as described in K. P 
aark.Exu-acting new information from the shape of the blood pressure pulse. Master's thesis, 
Massachusetts Institute of Technology, Cambridge. MA, 1991. 
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Figure 5 shows .he modified Windkessel model. The aona and major arteries axe modeled as a 
sin^e elasuc chamber (« which siores .he blood ejeced from UK left ve„.ricle durmg a s^lol^ 
Th! distal vessels a« modeled as capaci.ive (C,) and resisrive («,) elements ,hron^ whtch the , 
blood drains during a diastole. The oscma.or, effect of blood propagation is taken m.o accoum 
5 by introducing an effecUve mass (/,). The dynamic equation for the upstream .s derrved as 
below, where Qc is ihe cardiac output: 

dt // 

J 

where Qocan be solved from the constitutive law of the arterial wall on the 7" node of the local 
model derived in the previous section: 

p ^^U S -IL2iZ£s.^JsJI) (27) 



Downstream Blood Flow 

20 Similarly, the downstream dynamics extends the distal boundary (f^, .» ^ ^^^^^^ 
Veins <L be easily modeled as a reservoir when concemed with anenal hemodynMn.cs. Stnc. 
IX "cries are being moni.ored. which is close .o v«ns. .he ineru. .enn in *e downs«am .s 
negligible. The classic Windkessel model is used .o model Ae downstream as shown ,n F.g. 6 
wle O is the compliance of *e vessels in downsueam. «^ is .he characensuc reststamx, 

25 is peripheral resisrance. P. is an effort source if no interest in venous dynam.cs. 

The dynamic equaUon for the downstream can be written as: - ■ . . 

Where S^can be sdved f«»n Ae algebraic equaUon and Ox, constituUve law of the arterid wrfl 
on the TV* node: 



wo 00/72750 



PCTAJS00/1S006 



11 



Entire Arterial Model 

In this section, the models for the local arterial hemodynamics and the upstream/downstream 
5 dynamics, described above, are integrated to represent an entire systematic arterial stream. . 



The entire arterial model has (2N+3) slate variables and two inputs, as defined as following: 



10 



u = [Qc,PvY :(2xl) 



i>J^ :(2N + 3)xl 



(30) 
(31) 



From the continuity equation given by (16) and the constitutive law of the arterial wall given by 
(20), the pressures P, can be expressed in terms of the above state variables as: 

15 For further analysis of the nature of the hemodynamic behavior of the arterial flow, we linearized 
the dynamic model for local arterial segment given in (22) and (23) as follows: 



dt 



dt Az 



(34) 



20 From the dynamics equation for the upstream eqs. (24) - (27), the downstream cqs. (28) - (29) 
and the local arterial segment eqs. (32) - <34), a state-space representation of the extended model 
can be described in the following format: 



x = Ax-¥Bu 



(35) 



25 where A and B are: 



>V-3>H2N*y) 



:(2N + 3)x(2Ar + 3) 



(36) 



PCT/USOO/15006 

WO 00/72750 

12 




Design of Kalman Filter 

5 Kal-na. fiUe. ^ popularly .o cs,™a« unknown .ute variables .ha. ^ ™a 

above hemodynamic sys.en>. ,he ob«rva.ion equaUon n,us. be .^'^ 
i„suun,en,aUon me,hods .o be used. As s.a.cd previously, .he ob,ec„ve of *e ~ «^ »^ 
c»,Un«ously es.ima« *e blood pressu« n««ly fron, noninvasive '^^""^'^^.f ^^'""^ 
,« a periphe^l skin surface, to accordance wi.h p«fen«l embodunenu. . K^man fil^ « des.gn«J 
bL on an elecoical impedance pleAysmograph (EIP) and .wo pho.ople*ysmographs. 

A pho.op.eU,ysmograph employs a pair of LED and pho.ode.ec.or .o moni.« ^J'^^^f^ 
ar.eria) diame.er. Suppose .ha. a pho.opleAysmognph is anached on .he stan surface over each 
,3 ri bTends of .^Tanerial segmen. under consideraUon. "n-en. 0. .wo observauon funcuons 
y, and » can be simply described as a funcUon of Ume by using s..u vanables as: 

y,(i) = S,(»). y2«) = Si.<') 

20 EIP uses four elecuodes .o measure ,he elecuical imped»,ce of .he a.«rial segmen. surrounded 
^.riuodes. EIP is known .o provide O-e absolule measuremen. o, vo-umonc change of *e 
Irial segmen.. -n^refore. supposing *a. *e elec.«les are loca.ed a. U,e 
...erial secern under conside«uon. 0. ompu. of EIP y, can be descnbed m «rms of 0,e s.a« 

variables as: 



25 



^3(0 



Denning y<0 = b.VXy.(.').yM observaUon e,uaUon finally be defined as 



i 



I 



I 
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yO) = Cx{t) 



(38) 



where 
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Since a process noise and a measuremenl noise inherently exist, the state equations given in 
Section 2.4 must be extended as: 



10 



y = Cx + w 



(39) 
(40) 



where v and w are white noise processes, having known spectral density matrices, V and W, 
respectively. 

15 Using the above equations, the state variables x(t) can be estimated by the following dynamic 
equations: 



i = Ajc + Bi4 + /f(y-y) 



(41) 
(42) 



20 



where >>(/) is the estimated measurement, i(/) is the estimated state variables, and K is the 
Kalman gain matrix, which is updated as: 



25 



(43) 
(44) 



where M (r) is the covariance matrix of the state estimation error x(t) = x{t) - 5(r) . In the above 
derivation, we assume that v and vv are uncorrelated. By updating the Kalman gain based cm the 
nature of the process noises as described in the above equation, the Kalman filter provides the 
30 optimal estimation of the state variables. Finally, the internal blood pressures PfO) can be 

estimated by substituting the estimated state variables into (32) as: 
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p^=j£p^(S,-X(e,-e,.,)-So) /or,=1.2.-,Ar (45) 

The main issue in designing the above Kalman filter is whether the system given in (39) and (40) 
is observable or not. If the system is not observable, a Kalman filter can not be constructed to 
5 estimate the whole state variables. As it is found in the next section, the above system is not 
observable. However, the observability analysis to be provided in the next section will prove that 
the blood pressure given in (45) can be estimated from an observable subspace of the system. 



10 



Observability Analysis 



Observability Test 



t 

There are many criteria for testing the observability of a system, as described in W. S. Spector, 
"Handbook of Biological Data". Philadelphia Publisher. 1956. The standard test is the 
15 "Algebraic Controllability Theorem," as described by T. Kailaih, Unear Systems, Prentice-Hall, 
NJ, 1980, and it simply states: 

A system {A.C) of order n is observable if and only if the rank of the observability test 
matrix 

o=lc^i4'"c^.••.(A'^)-c^] (46) 

is equal to n. 



20 



This is arguably the easiest criterion to test the observability of a system. 

25 The above observability test was applied to the (2N+3)-th order system given by (35) and 
(38), and it was found that the rank of the observability matrix is 4 when N=3 or 3 when N>3, 
which is smaller than the order of the system. Therefore, the system is not observable and a state 
estimator such as a Kalman filter cannot re-construct the whole state variables. However, it wiU 
be found that the blood pressure given in (45) can be estimated from a part of the state variables 

30 and the part-lies in the observable subspace. of the state space. Namely, the blood pressure c^ be 
estimated from a set of the state variables which are observable with the Kalman filter designed 
in the previous section. To prove this argument, the whole state variables are decomposed into an 
observable sub-space and an unobs rvable sub-space. 
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ObservableAJnobservable Sub-space Decomposition 

A staircase algorithm is used for the state-space decomposition. Letting r to be the rank of 
the observability matrix given in (46), for the system described by (35) and (38), there exists a 
5 de-coupling similarity transformation matrix T such that 



A = TAT'^ = 



no 



+ 3-r)>.(2W + 3-r) A,, : <2*f * 3 - r) » H 
0 A„:r*t J 



C=Cr=l0 C„:3x,] 

Namely, the state equation can decomposed into an observable subspace and an unobservable 
subspace, and the r-dimensional observable subspace is represented by IA„ , Co]. Suppose T is 
expressed as 

10 I 

7' = [7,/ :«W*3-Oxt2« + 3) r/ :rx«Ar + 3)J ^^^^ 

Then, the transformed state variables z are decomposed into the observable state variables and 
the unobservable state variables Zuo as: 

r T If I 

(49) 



=1 z„='- J" ki 



Consequently, the set of the transfoimcd state variables Zc=ToX is observable from the output 
given in (38) using the Kalman filter designed in the previous secuon. 

20 

Blood Pressure Estimation from Observable Sub-space 

The blood pressure can be calculated from state variables, according to eq. (45). This 
equation can be expressed in a vector form such as: 

where G, is a (2N + 3) x 1 row vector. 

From the state-space analysis, it is found that there exists a rxl row vector H such that 
G=HTo. Therefore, the Wood pressure in (50)ean be described as 

P,^G,x = HT,Jc = Hz„ (51) 



* 
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Namely, the blood pressure can be estimated from the observable variables Zo- 

The above analysis of the state-space decomposition shows thai two Photoplethysmographic . 
sensors and one EIP sensor on an arterial segment can estimate the pressure wavefonns using a 
5 Kalman filter. Based on these results, a cuff-less ambulatory blood pressure monitoring device 
can be designed. 

Figure 7 illustrates the sensor configuration. Photoplethysmographic sensors 10 and EIP 
sensors 20 are located on the two bands, which fit a human finger. The telemetry 30 can transmit 
signals wirelessly. 



10 • 



Simulation 



Numerical simulations have been conducted to verify the approach. The hemodynamic 
process was simulated using MATLAB on a PC, and the Kalman filter for the blood pressure 
15 estimation was applied to the simulated process. The pressure estimated by the Kalman Filter 
was compared with the digital blood pressure measured by an arterial tonometer. 

Simulation Setup 

20 The simulation was conducted using hemodynamic parameters of a digital artery because 
many finger plethysmographs are commercially available and easy to be miniaturized. The 
following parameter values were used for the simulation: 

Blood density p = L06 gr/cm^ , 
25 Blood viscosity fi = 0.04 poise , 

Radius of digital anery r = 0.5 mm , 

Arterial wall viscosity 77 = 700 dyn-s/cm^ , 

Arterial wall elastic modulus E^7xl& N/m^ , 

Characteristic resistance Red = Llxl(f dyn-s/cm^ ^ 
30 Peripheral resistance Rpd = 1.2xl(f dyn s/cm^ , 

Downstream compliance Cd = LlxlOT^ cm^ /dyn , 

Length of digital artery segment L = 7 cm. 

Nodes of the system N = 3. 

35 The above parameter values were obtained from published literatures such as E. Belardinelli 
and S. Cavalcanti, "A New Nonlinear Two-Dimensional Model of Blood Motion in Tapered and 
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Elastic Vessels," Comul. Biol. Med., Vol. 21, pp. 1-13, 1991; W. S. Spector, "Handbook of 
Biological Data", Philadelphia Publisher, 1956; B. M. Leslie, et el., "Digital Artery diameters: 
An anatomic and clinical study". Journal of Hand Surgery. Vol. 12A, No. 5, Part 1, pp740-743, 
Sep. 1987; H. Power, "Bio-fluid Mechanics", Computational Mechanics Publications, Boston, 
5 1995; and KJ. U, "Arterial System Dynamics", New York University Press, New York, 1987. 
The distributed model of a digital artery used in simulation is shown in Figure 8. For the 
simplicity of the simulation, the upstream dynamic in the arterial hemodynamic model was not 
included. Instead, measured blood pressure signals at Section Si were used as an input to the 
serial model. The definitions of the inputs, state variables and outputs in this simplified model 
10 are 

inputs [Pj Pv] • 

7* 

Slate variables x = [Qj Q2 Sj S2 Ss Pdl » 
outputs y = /Si V Ss] . 

15 In this setup of simulation, similarity transformation matrix J can be calculated 

numerically: 
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where the last four rows of T matrix represent the observable subspace To. 

20 

Input Pi is measured by an arterial tonometer (MILLAR, TX). The other input, venous 
pressure Pv is assumed as a constant (20mmHg). A profile of input Pi is shown in Figure 9. 

Outputs Sj and S3 are measured by a pulse plethysmograph (CB Sciences, Dover, NH) 
25 and Vis measured by an electrical impedance plethysmograph (Parks Medical Electronics, 
Aloha, OG). The measurements are shown in Figure 10. 



Simulation Results 



tit 
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The Kalman filler consiructed in Section 3 is simulated in MATLAB to estimate slate 
variables and blood pressure. Pj, Pv (u in eq. (41), depicted in Fig. 9) and measurement Yj, Y2, Y3 
(y in eq. (41), depicted in Fig. 10) are feed into a Kalman Filter. The Error covariance and' 
Kalman filter gain art calculated for each sample of the sequence and state variables are updated 
5 according to eq. (43) and (44). Necessary state variables are then substituted into eq. (51) to 
estimate blood pressure. 

Figure 1 1 shows the comparison between the measurement and Kalman Filter estimation 
of the output, in which it can be seen that the Kalman Filler works very well to reduce white 
10 Gaussian noises as expected. 

Figure 12 shows the comparison between the measurement and the Kalman filter 
estimation of blood pressure. 

15 From the results shown in Fig. 1 1 and Fig. 12, it can be concluded that a Kalman filler is 

very robust to noise, especially white noise. It is feasible to estimate blood pressure accurately 
based on the measurements from plethysmographs and a hemodynamic model. 

The Monitoring System 

20 

An exemplary monitoring system is shown in Figure 13. Outputs Si and ^3 are measured 
on the left hand middle finger by dual photo plethysmograms JO and Vis measured on the same 
finger by an electrical impedance plethysmogram 20. A conu-oller 40, not necessarily consisting 
of any of the elements shov^m, derives a measure of the blood pressure. Prior to operating, the 
25 system is calibrated against a pressure cuff or other blood pressure monitoring device. 

Although various exemplary embodiments of the invention have been disclosed, it should 
be apparent to those skilled in the art that various changes and modifications can be made which 
will achieve some of the advantages of the invention without departing from the mie scope of the 
30 invention. These and other obvious modifications are intended to be covered by the claims that 
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What is claimed is: 

1. A moniioring sysiem for moniioring a blood pressure of a patient, the monitoring 
system comprising: 

a. a first phoiopleihysmograph proximate to a finger of the patient for providing a 
measure of change in the anerial diameter at a first location of a specified anery 
of the patient; 

b. a second phoioplethysmograph proximate to the finger of the patient and 
displaced relative to the first phoioplethysmograph for providing a measure of 
change in the anerial diameter at a second location of the specified artery of the 
patient; 

c. an electrical impedance plethysmograph in electrical contact with the finger of 
the patient for providing a measure of change in the electrical impedance of an 
anerial segment between the first and the second locations of the specified 
anery; and 

d. a controller deriving a measure of the blood pressure of the patient based on the 
measures of change in the anerial diameter at the first and second locations of 
the specified anery and the measure of change in the electrical impedance of an 
anerial segment. 

2. A monitoring sysiem according to claim 1 , wherein the first phoioplethysmograph is 
bome by the patient on a finger ring. 

3. A moniioring sysiem according to claim 1, wherein the first phoioplethysmograph is 
bome by the patient on a first band of the finger ring and the second 
phoioplethysmograph is bome by the patient on a second band of the finger ring. 

4. A monitoring sysiem according to claim 1, funher including a transmitter for 
transmitting the measure of the blood pressure of the patient to a remote location. 

5. A monitoring sysiem for monitoring a blood pressure of a patient, the monitoring 
system comprising: 

a. a monitor having a first and a second band to be worn by the patient on a single 
finger, the monitor comprising: 

i. a first phoioplethysmograph disposed on the first band for 

providing a first signal based on a first anerial diameter of the 
patient; 
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ii. a second phoiopleihy sinograph disposed on the second band 
for providing a second signal based on a second arterial 
diameter of the patient; 

iii. an electrical impedance plethysmograph disposed on the first 
and second bands for providing a third signal based on the' 
electrical impedance of the a segment of an artery of the 
patient; 

and 

b. a controller for analyzing the first, second , and third signals and detcnnining a 
measure of the blood pressure of the patient. 

6. A method for monitoring the blood pressure of a patient, the method comprising: 

a. deriving a measure of change in the diameter of a first end of a segment of an 
artery of the patient; 

b. deriving a measure of change in the diameter of a second end of the segment of 
the anery of the patient; 

c. deriving a measure of the volume of the segment of the anery of the patient; 

d. applying a model of arterial blood How to the derived measures of change in the 
diameters of the first and second ends of the arterial segment and the volume of 
the segment for calculating the instantaneous blood pressure of the patient. 

7. A method in accordance with claim 6, wherein the step of deriving the measure of 
change in the diameter of the first end of a segment of an artery includes receiving a 
signal of a first photoplethysmograph. 

8. A method in accordance with claim 6, wherein the step of deriving the measure of 
change in the diameter of the second end of a segment of an anery includes 
receiving a signal of a second photoplethysmograph. 

9. A method in accordance with claim 6, wherein the step of applying a model includes 
applying a Kaman filter for estimating internal state variables. 
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